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Abstract 

We examine the lepton flavor-changing processes in i?-parity violating MSSM. First, we update 
the constraints on the relevant i?-violating couplings by using the latest data on the rare decays 
ii — > ij'y. We find that the updated constraints are much stronger than the old ones from rare 
Z-decays at LEP. Then we calculate the processes Z — > iiij and 77 — > £{£j . We find that with the 
updated constraints the i?-violating couplings can still enhance the rates of these processes to the 
sensitivity of GigaZ and photon-photon collision options of the ILC. 
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I. INTRODUCTION 



The minimal supersymmetric model (MSSM) is a popular extension of the Standard 
Model (SM). In this model the invariance of -R-parity, defined by R = ^—±Y+ 3B+L f or a field 
with spin S, baryon-number B and lepton-number L, is often imposed on the Lagrangian in 
order to maintain the separate conservation of baryon-number and lepton-number. Although 
i?-parity plays a crucial role in the phenomenology of the MSSM (e.g., forbid proton decay 
and ensure a perfect candidate for cosmic dark matter), it is, however, not dictated by 
any fundamental principle such as gauge invariance and there is no compelling theoretical 
motivation for it. The most general superpotential of the MSSM consistent with the SM 
gauge symmetry and supersymmetry contains R- violating interactions which are given by [l| 
1. _ 1 



W# - -XijkLiLjEl + M ijk LiQjD c k + -\" jk e abd U- a D c jb D c kd + /^Ljif 2 , 



(1) 



where k are generation indices, c denotes charge conjugation, a, b and d are the color 
indices with e abd being the total antisymmetric tensor, H 2 is the Higgs-doublet chiral super- 
field, and Li(Qi) and Ei(Ui, Di) are the left-handed lepton (quark) doublet and right-handed 
lepton (quark) singlet chiral superfields. The dimensionless coefficients A^ (antisymmetric 
in i and j) and \' i j k in the superpotential are L- violating couplings, while A" jfc (antisym- 
metric in j and k) are I?-violating couplings. So far both theorists and experimentalists 
have intensively studied the phenomenology of i?-parity breaking supersymmetry in various 
processes [2J, |3j and obtained some bounds [4|. 

The lepton flavor- changing (LFC) processes, which have been searched in various exper- 
iments 5H3], are a sensitive probe for new physics because they are extremely suppressed 
in the SM but can be greatly enhanced in new physics models like supersymmetry In 
-R-parity breaking supersymmetry, these rare processes may receive exceedingly large en- 
hancement since both A and A' couplings can make contributions. Such enhancement was 



considered in the decays k — > ^7 [9( and Z — > £i£j [1 01] , the \i — e conversion in nuclei ll| , 
and the di-lepton productions pp/pp — > e ± /i =F + X [12] and e + e~ — > e ± /i =F 131 ] . 

Since the GigaZ and photon-photon collision options of the ILC can precisely measure 
the LFC processes Z — > £i£j and 77 —> £i£j (i 7^ j and £i = e, /x, r), we in this work study 
these processes in i?-violating MSSM. Noting that the experimental upper bounds on the 
LFC r-decays became more stringent recently 6], we will first check the constraints on the 



relevant i?-violating couplings from the latest measurement of 



j7. Then, with the 



updated bounds on the relevant i?-violating couplings, we calculate Z — > £i£j and 77 —> £i£j 
to figure out if they can reach the sensitivity of the GigaZ and photon-photon collision 
options of the ILC. 

The paper is organized as follows. In Sec. II we describe the calculations for £^ — > £jj, 
Z — > £i£j and 77 — > £i£j. In Sec. Ill we present some numerical results and discussions. 
Finally, a conclusion is drawn in Sec. IV. 



II. CALCULATIONS 

In terms of the four-component Dirac notation, the Lagrangian of the L-violating inter- 
action is given by (in our calculations we take the presence of \'^ k as an example) 



C\> — —X'ij k 



V 



L44 + 444 + (4)*(4) c 4 



-zi44-«14^-(4mM +h.c. (2) 



The LFC interactions £{£jV (V = 7, Z) are induced at loop level by exchanging a squark u 3 L 
or d k R , which is shown in Fig.l. 
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FIG. 1: Feynman diagrams for — £a transition induced by the L- violating couplings at one-loop 

level. 

For the decays £i — > £j r y we take // — > as an example to show the analytic results. The 
gauge invariant amplitude of \i — > e7 is given by 

M(/i -> e 7 ) = 2Aw(p e )P Ji (2e-p A( - 4)u{p„), (3) 
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FIG. 2: Feynman diagrams for 77 — > l^lj induced by the L-violating couplings at one-loop level. 
The effective 7 — l\ — £j vertex in (a,b) is defined in Fig. 1. 



where A is given by (assuming the degeneracy for squark masses) 
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The decay branching ratio reads 

BR(fj, -> e7) 



48vr 



(4) 

(5) 
(6) 

(7) 



For the decays Z -> ^ we calculate the decay rates numerically by using the effective 
vertex presented in Appendix A. Note that according to the effective vertex method [14J, 
the external legs of the effective vertex can be on-shell or off-shell and thus the vertex can be 
used in any relevant process. The expression in Eqs.(3-6) can be obtained from the effective 
vertex in Appendix A by putting both leptons on shell. 

For the process 77 — > £{£j, besides Fig. 2 (a,b) induced by the effective vertex given 
in Appendix A, more diagrams shown in Fig.2 (c-i) also come into play. The analytic 



expressions of the amplitudes of these diagrams are given in Appendix B. These amplitudes 
contain the Passarino-Veltman one-loop functions, which are calculated by using LoopTools 



261 ] . We checked that the amplitudes have gauge invariance and the ultraviolet divergence 



■^77->£i£j — / dy/s^ - T!_cr 11 ^£.£.{s 11 ) = £ e + e -cr 77 ^.^(s) (8) 



cancelled. 

Since the photon beams in 77 collision are generated by the backward Compton scattering 
of the incident electron- and the laser-beam, the events number is obtained by convoluting 
the cross section of 77 collision with the photon beam luminosity distribution: 

d£>- 

s 77 

where dC^/dy/s is the photon-beam luminosity distribution and cr^^^i (s) ( s is the 
squared center-of-mass energy of e + e~ collision) is defined as the effective cross section of 
77 — > iilj. In optimum case, it can be written as [15[ 

Vry-Hiiji 3 ) = / 2zdza^(i.z.(sy 7 = z 2 s) — F 7/e (x)F 7/e (— ) (9) 

J y/a J z ' 2 / X max X X 

where F 7 / e denotes the energy spectrum of the back-scattered photon for the unpolarized 
initial electron and laser photon beams given by 



F l/e (x) 



D(0 
with 



1 Ax Ax 

1-X+- — r + 



l-x £(l-x) i 2 {[l-xf 



(10) 



Here ^ = AE e E Q /m1 (E e is the incident electron energy and E is the initial laser photon 
energy) and x = E/Eq with E being the energy of the scattered photon moving along the 
initial electron direction. 



III. NUMERICAL RESULTS AND DISCUSSIONS 



In our calculations we take the SM parameters as 16] 



m M = 0.106 GeV, m T = 1.777 GeV, m b = 4.2 GeV, a = 1/137, sin 2 9 W = 0.223 (12) 



The top quark mass is taken as the new CDF value mt = 172.3 GeV [17]. The relevant 
SUSY parameters in our calculations are the masses of squarks as well as the i?-parity vio- 
lating couplings listed in Table I. The strongest bound on squark mass is from the Tevatron 



experiment. For example, from the search of the inclusive production of squark and gluino 
in i?-conserving minimal supergravity model with Aq = 0, fi < and tan/3 = 5, the CDF 



gives a bound of 392 GeV at the 95 % C.L. for degenerate gluinos and squarks |18|] . However, 
this bound may be not applicable to the i?-violating scenario because the SUSY signal in 
case of R- violation is very different from the i?-conserving case. The most robust bounds on 
sparticle masses come from the LEP results, which give a bound of about 100 GeV on squark 

n 

or slepton mass |19j. In our numerical calculations, we assume the presence of the minimal 
number of R- violating couplings, i.e., for each process only the two relevant couplings (not 
summed over the family indices) are assumed to be present. 
For £i — > £j7, the latest experimental data is j3] 

BR(fj, -»■ &y) < 1.2 x 10~ n , (13) 
BR(t &y) < 1.1 x HT 7 , (14) 
BR(t -> pi) < 4.5 x HT 8 . (15) 

We use these data to update the bounds on the relevant L-violating couplings. The new 
bounds are compared with the old ones in Table I for rriq = 100 GeV (here we take squark 
mass of 100 GeV for illustration and for heavier squarks the bounds on the L-violating 
couplings will become weak, as will be shown later). We can see that the new bounds are 
much stronger than the old ones. Since the bounds on X' i33 Xj 33 (i ^ j) are weakest, we only 
consider the contribution of A^gA'-gg in our following numerical calculations. 



Note that the neutrino masses could also constrain the A' couplings, especially A^ 20]. 
But these constraints depend on more parameters in addition to the squark mass. For 
example, the one-loop A' contributions to the neutrino masses are sensitive to the left-right 
squark mixings and the two-loop contributions further involve the slepton mass. For small 
squark mixings with appropriate sign, there may exist a strong cancellation between one- 
loop and two-loop effects, and in this case, the constraints from the neutrino masses can be 
avoided. Since the aim of our study is the sensitivity of the LFC processes to A' couplings 
and the A' contributions to these LFC processes are irrelevant to the additional parameters 
involved in the contributions to the neutrino masses, in our analysis we did not consider 
such constraints from the neutrino masses. 
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TABLE I: Our new upper bounds on the 
data jjj, in comparison with the old ones 4} 



L-violating couplings for nig = 100 GeV from 



couplings 



\' X' X' X' 

A 111 A 211> A 112 A 212 

y x' 

A 113 A 213 

x' v v x' 

A 121 A 221> A 122 A 222 

X' X' 

A 123 A 223 



X' X' 

A 131 A 231' 

\' X' 

A 133 A 233 



X' X' 

A 132 A 232 



X' X' X' X' 
A 111 A 311> A 112 A 312 

X' X' 

A 113 A 313 



A 121 A 321> 

X' X' 

A 123 A 323 

X' X' 

A 131 A 331> 

X' X' 

A 133 A 333 



A 122 A 322 



y X' 
A 132 A 332 



X' X' X' X' 
A 211 A 311' A 212 A 312 



X' X' 

A 213 A 313 

X' X' 

A 221 A 321> 

X' X' 

A 223 A 323 

A 231 A 331> 

X' X' 

A 233 A 333 



X' X' 

A 222 A 322 



A 232 A 332 



New bounds Old bounds [4] 

5.7 x 10~ 4 
5.7 x 10~ 4 



7.74 
7.85 
7.78 
7.89 
1.27 
1.63 
5.54 
5.56 
5.57 
5.65 
9.06 



1.17 x 10~ 2 

3.55 x 10~ 4 

3.60 x 10~ 4 

3.56 x 10~ 4 

3.61 x 10~ 4 
5.80 x 10~ 3 
7.48 x 10~ 3 



x 10~ 5 
x 10~ 5 
x 10~ 5 
x 10~ 5 
x 10~ 3 
x 10~ 3 
x 10- 4 
x 10- 4 
x 10- 4 
x 10- 4 
x 10~ 3 



5.7 x 10~ 4 
5.7 x 10~ 4 
7.7 x 10~ 3 
1.0 x 10~ 2 
1.2 x 10~ 2 
1.2 x 10~ 2 
1.2 x 10~ 2 
1.2 x 10~ 2 
1.2 x 10~ 2 
1.2 x 10~ 2 



For Z — > £{£j, the upper limits from LEP are 2jJ, |22] 



BR(Z lie) < 1.7 x 1(T 6 , 
BR(Z re) < 9.8 x 1(T 6 , 
BR(Z ->■ th) < 1.2 x 10" 5 . 



(16) 
(17) 
(18) 



The bounds from these LEP data are compared with the bounds from £i — > £jj in Fig.3. 
One can see that the upper bounds on the couplings from the LEP Z-decay data 2l|, |22|] are 
weaker than the ones from L — > £^ data [7] . Note that the bounds from the LEP Z-decay 
data were also studied in [lOj and our results are consistent with theirs except that in 10] 
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FIG. 3: Various bounds on the L-violating couplings versus the squark mass. The solid, dashed 
and dotted curves are the bounds on A'^A^, A133A333 and A233A333, respectively. Also shown are 
the 2<T sensitivity from Z-decays at GigaZ and the 3<t sensitivity from 77 — > e(or fi) r at the ILC 
with center-of-mass energy of 500 GeV and a luminosity of 3.45 x 10 2 /& _1 . 



the sum over index k is implied for A'^A^ with rriq = 200 GeV. 

The possible sensitivity of GigaZ to the LFC decays of Z-boson could reach 

BR(Z -> fie) ~ 2.0 x 10~ 9 , 
BR(Z -> re) ~ k x 6.5 x 10~ 8 , 
BR(Z -> Tfjt) ~ k x 2.2 x 10~ 8 



23] 



(19) 
(20) 
(21) 
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with the factor k ranging from 0.2 to 1.0. In Fig. 3 we take k = 1.0 to show the sensitivity. 
In contrast to the i?-conserving case in which only Z — > fir is accessible at the GigaZ jg], 
the i?-violating couplings under the bound from /j — > can still enhance all the channels 
Z — > £i£j to the sensitivity of the GigaZ. This implies that the GigaZ can further strengthen 
the bounds on A^ 33 A^ 33 in case of un-observation. These bounds, unlike the constraints from 
neutrino masses which involve more parameters, are only dependent on the squark mass. 
For the 77 collision results shown in Fig. 3, we fixed the parameters as £ = 4.8, D(^) = 



1.83 and x max = 0.83 [15|. Since the L- violating couplings relevant to the process 77 — > ep 
is stringently constrained by // — > ej, we in Fig. 3 only show the results for the channels 
with a tau lepton in the final states, i.e., 77 —> ef, /if. The background for 77 — > ef 
comes from 77 — > t + t~ — > T~y P p T e + , 77 — > W + W~ — > T~v e v T e + and 77 — > e + e~r + r~, 
and we make kinematical cuts [l3j: | cos6^| < 0.9 and p l T > 20 GeV [I = e,/i), to enhance 
the ratio of signal to background. With these cuts, the background cross sections from 
77 — > t + t~ — > T~v e v T e + , 77 — > W + W~ — > T~v e v T e + and 77 — > e + e~r + r~ at y 7 ? = 500 
GeV are suppressed respectively to 9.7 x 10 -4 fb, 1.0 x 10 _1 fb and 2.4 x 10~ 2 fb (see Table I 
of [3]). To get the 3cr observing sensitivity with 3.45 x 10 2 fb -1 integrated luminosity 24], 
the production rates of 77 — > ef,fif after the cuts must be larger than 2.5 x 10~ 2 fb 13]. 
We see from Fig. 3 that under the current bounds from li — > ^7, the L-violating couplings 
can still be large enough to enhance the productions 77 — > ef, fif to the 3a sensitivity. 

We also show the cross sections of 77 — > iilj as a function of center-of-mass energy ^fs 
of the ILC in Fig. 4. We see that with the increasing of the center-of-mass energy, the cross 
sections of these processes become smaller. Such a behavior is similar to the results in the 



R-conserving MSSM shown in 13]. 



Finally, we point out that the LFC processes can also put bounds on the products A^ 31 Aj 31 
and A / i 3 2 A^3 2 , and our numerical results indicate that such bounds are quite similar to those 
in Fig. 3. We note that these bounds on A£ 31 A'- 31 and A / i 32 A'- 32 from Z — > klj at GigaZ are 



generally stronger than those from the neutrino masses 20j 



IV. CONCLUSION 

We evaluated the lepton flavor- changing processes in i?-parity violating MSSM. First, we 
used the latest data on the rare decays li — > ijj to update the constraints on the relevant 
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FIG. 4: The cross sections of 77 — > as a function of center-of-mass energy yfs. The couplings 
A / 1 3 3 A / 233, A' 1 33A 3 33 and A233A333 are fixed at their upper bounds at Mg = 300 GeV. 

R- violating couplings. Then we calculated the processes Z — > 1$ j and 77 — > £i£j. We found 
that with the updated constraints the i?-violating couplings can still enhance the rates of 
these processes to the sensitivity of GigaZ and photon-photon collision options of the ILC. 
So, the GigaZ and photon-photon collision of the ILC can either observe these A'-induced 
LFC processes or further strengthen the bounds on the A' couplings in case of un-observation. 
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Appendix A: Expressions of effective vertex j(Z) — l% — £j 



Here we list the expression for the L- violating contributions to the effective vertex j(Z) — 
e — fx. Other effective vertices "f(Z) — 11 — r and ^y(Z) — e — r are similar to j(Z) — e — li 
and can be obtained by replacing the corresponding momentum and mass. The the effective 
vertex j(Z) — e — li is given by 

r j{Z)e^ = r l(Z)e^~L} + if*)** , (Al) 

where the two terms on the right side denote the L-violating loop contributions by exchang- 
ing respectively the squarks and d^, given by 

rTC^Pe)^ = ae{~[Cl pl a lx^ ~ Cl^,-^ e ) lxl a ]PL + \[2C 2 aXl a - Cl{p,+ Pe ) xl a ]P L 

-X^-fB^PL + \\l a i,lxPL + m^ a lx P R )]Bl - \m 2 dk Ch x P L \ (A2) 

rr(P,.P e )b-« = ™{-\{Clpl a lxl P - Cl(^-i e ) lxl a )P L + \[2C&f - Ct(p, +Pe)xl a ]P L 

-\lH»l a BlPL + \\l a ielxPL + m^ a lx P R )}Bt - \m 2 Ch x P L \ (A3) 
m fi m ti 6 ) 

r Z x e ^,p e )U f = 6e|^[^ /37 a 7A7 /3 -^(^-Ve)7A7 a ] J PL-(l-^)mi fc C 5 7Aj PL 

4s 2 1 - 2s 2 
+ (1 - -^)[2C 6 aXl a - Clip, +p e ) xl a ]P L ^jrfBlPL 

- 24) 7 ^ e7A Pi - 2s 2 w m^ a 7xPR)}B 6 a ) (A4) 
T Z x e ^ PlM p e )\ S n = 6e|-(l-^)[^ 7 W-^(^-^)7A7 a ] J PL + ^<.C 7 7APL 

-^c s aX r - ci(p,+p e ) x r]p L - ^-^-rtt-fBiPL 

+4^(1 - 2s 2 w hy e ^ x P L - 2slm^ a jxP R )]B 8 a ) (A5) 
m 2 J 

with a = ~3g~5^ ) b = ^t^^T anc ^ P e anc ^ 2V denoting respectively the momenta of 
the electron and muon. In the above expressions, the functions B % a and C % a a/3 are the 
Passarino-Veltman functions. For these loop functions, we adopt the definition in {25 1 and 



use Loop Tools [26j in the calculations. The functional dependence of these loop functions is 



If 



given by 



c\- 


2 2 2 \ ri2f 

-p^p e ,m dk ,m-L,m dk ), C (- 


-Pe,Pe 


-Pp.; 




,m 2 - L . 


u j 


(A6) 


c\- 


-P»,Pe,rnl.,m d : R ,ml j ), C 4 (- 


-Pe,Pe 


-Pix, 


3 


d k 


m %) , 

d k 


(A7) 


C 5 (- 


-p^,p e ,m 2 dk ,m 2 L,m 2 dk ), C 6 (- 


-Pe,Pe 


-Pli; 




, m 2 L . 




(A8) 


C 7 {- 


-pu,Pe,ml.,m 2 , R ,ml ), C 8 (- 


-Pe,Pe 


-P/J., 


m u 3 ; 


mi, 


d k ' 


(A9) 


B\ 


-Pn,rn 2 dk ,m 2 ^), B 2 (-p e ,m 2 dk . 












(A10) 


5 3 (- 


-P^rn: I ,m} R ), B 4 (-p e ,m 2 

k J 




) 








(All) 


5 5 (- 


-p„m 2 dk ,m 2 ,), B\-p^m 2 dk 












(A12) 


5 7 (- 


-P^m 2 Uj ,m 2 iR ), B 8 {-p e ,m 2 Uj , 


,m] R ) 










(A13) 



Appendix B: Expressions of amplitudes for 77 — > £i£j 

The amplitudes of the diagrams in Fig.2(a-i) are given by 

M (a)\uj,d* = ^e)(iej x ) ^^ r], e ^(-p^p2 - p^L dR v(fi)e^e p 2 (Bl) 

M (b)\uf,d* = u{e)Tl e ^{ P2 - Pt,,Pe)\ U L :d R ^ % _ j-^ {ie~i p )v{ii)e\e p 2 (B2) 
i 8 11 

M(c) '"i = ~Y^2^Q e2 ^ X ^k X '3jku{e)Cl-f a P L v(fi)eve 2 (B3) 
M(c)l ^ = I^(r 2 ) A ^4^( e )^°^ P ^^) e i- e 2 (B4) 
M (d) | fi i. = I ^(^e 2 )A' ljfc A; jfc n( e )|^ /35 7 Q 7p7 /3 7A7 <5 + ^ / 37 Q 7p^27A7 /3 

+^(^1 +y 2 )7p7 Q 7A7 /3 + +^2)7p^27A7° + "^^TpTA 

+7p7A7° + 7p7°7a) + m^Z^i +^ 2 )7p7a + 7p^27A]}i 3 L^)ei4 (B5) 

-2DL(2p e -Pi)p + ^(2Pe-Pi)p(2pe " 2 Pl - p 2 ) j 7 a P L t;^)ey^ (B6) 
= ^(-^e^A'^A^e)!^^ _ D%{2p e - p 2 ) x ] 1 a lp ^ - [2D 3 aX 

-D 3 a (2 Pe ~ P2)a]7°7p^i + m 2 dk [2D\ - Dl(2p e - p 2 ) x ] 7p ) P L v(p)e^ 2 (B7) 
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M i 9 M = i^(^ e2 )4fe4fc«( e ){^7 Q 7p7 /3 7A7' 5 + ^7 a 7p^7A7 /9 

+D£pWi +^ 2 )7p7 a 7A7 /? + £>a(Vl +^2)7p^27A7 Q + ^^( 7 a 7p7A 
+7p7A7° + 7p7"7a) + ^i^oK^i +^) 7 p7A + 7pMa] ^P L v{p)e\e p 2 (B8) 

= Y^i-^ 2 )^ k ^ jk ^e)UDl Xa -2Dl a (2p e -2 Pl -p 2 ) x 

-2D\ a {2p e - Pl ) p + l£(2p c -pi) P (2 Pe - 2p x -p 2 )\ 1 a P L v{ l i)e\e p 2 (B9) 

-l£(2p e -P2)a]7 q 7p^i +m%[2Di - D§(2p e - p 2 ) x } 1( ^P L v(M4- (BIO) 
Here the effective vertices appearing in Eqs. flBll) and (1B2I) are defined in Appendix A. 
The amplitudes for the diagrams with the two photons exchanged are not presented here, 
which can be obtained from the above corresponding amplitudes with replacement p\ -H- p 2 
and e\ <-¥ e 2 . The functional dependence of the Passarino-Veltman loop functions C l a and 

C 9 (- Pe ,Pn + p e ,m 2 d ,m 2 ~ L ,m 2 ~ L ), C 10 (-p e ,p M + p e , ml ,mL ml) (Bll) 
D 1 (P2,Pi,~Pe, m\ , , m\ , m 2 L ) , L> 2 ( -p e , Pl , p 2 , m\ , m 2 L , m 2 L , m 2 L ) (B 1 2) 

3 3 3 3 

D 3 (-p e ,p 2 , -pp , m\ , m^L , m\t , ) , D 4 (p 2 ,px,-p e , m\. , , ml. , m 2 u . , m|« ) (B 13) 
D 5 (- Pe ,p 1 ,p 2 , m 2 u , , m\ R , m\ R , m\ R ) , D 6 (-p e ,p 2 , -p^ , m^. , m |h , m| fl , m 2 u . ) (B14) 
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